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In order to understand the origin of one-body dissipation in nuclei, we analyze the behavior of a 
gas of classical particles moving in a two-dimensional cavity with nuclear dimensions. This "nuclear" 
billiard has multipole-deformed walls which undergo periodic shape oscillations. We demonstrate 
oo : that a single particle Hamiltonian containing coupling terms between the particles' motion and 

the collective coordinate induces a chaotic dynamics for any multipolarity, independently on the 
geometry of the billiard. If the coupling terms are switched off the "wall formula" predictions 
are recovered. We discuss the dissipative behavior of the wall motion and its relation with the 
order-to-chaos transition in the dynamics of the microscopic degrees of freedom. 
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| In the last twenty years the dissipation of collective motion in nuclei has been widely observed |jj in low energy 
particle and heavy ion collisions and it still represents a theoretically unsolved problem. It is commonly believed that 
both one-body processes, i.e. collisions of nuclcons with the nuclear wall generated by the common self-consistent 
t-H , mean field, and two-body collisions produce dissipation, although their interplay is not well known 

The theory of damping and the approach to equilibrium of a collective (usually slow) degree of freedom coupled 
to non-collective (usually fast) degrees of freedom is a quite general one, both at classical and quantal levels. In the 
last decades several papers have been published on the subject, see for example refs. @-|Jll|,|l^,|l|-0. The most 
studied case is the extreme adiabatic limit, namely when the ratio between the characteristic times of the fast degrees 
of freedom and of the slow one is considered arbitrary small. Under the further assumption that the fast degrees of 
freedom have an ergodic dynamics for a fixed value of the collective variable, a description of the damping and approach 
to quasi-equilibrium can be given in terms of generalized diffusion equations . The theoretical framework for the 
treatment of this limiting case, at least at the classical level, seems to be well established, despite some questions 
appear to be still open Q. Unfortunately, in the majority of realistic cases the adiabatic condition is only partly 
satisfied and the ergodic condition is not often fulfilled. The damping, however, is still necessarily connected to 
some degree of chaoticity of the overall dynamics of the system. The connection between chaotic dynamics and 
dissipation in this more general situation has been studied by different authors but still lacks of a general theoretical 
framework. In particular Wilkinson M] presented some evidences which indicate that the occurrence of an integrable 



ON 



X 



or nearly integrable motion of the fast degrees of freedom strongly suppresses both the fluctuations and the speed of 
relaxation towards equilibrium. Related studies have been developed in ref. fTfl| , where a coexistence of chaoticity 
and slow relaxation to equilibrium in a Hamiltonian Mean Field (HMF) model has been found, and in ref. Jlq ), where 
the relation between chaoticity and the approach to equilibrium in a hard sphere gas and a Lorentz gas has been 
analyzed. 

In this context, Blocki et al. analyzed the behavior of a gas of classical non-interacting particles enclosed in a 
multipole-deformed container which undergoes periodic shape oscillations. Particles move on linear trajectories and 
collide elastically against the walls. In this model wall and particles' motion is uncoupled and therefore the wall keeps 
oscillating at the same frequency pumping energy into the gas. For this system, the authors study the increase of the 
particles' kinetic energy as a function of time. They find that for octupole and higher modes the gas kinetic energy 
increases with time, in agreement with the "wall formula" predictions (lfj. They attribute the different behavior to the 
fact that for low multipolarity deformation the particles' motion is regular and corresponds to an integrable situation, 
whereas for higher multipolarities the shape irregularities lead to divergence between trajectories and therefore to 
chaotic motion. Although their results look very interesting, their application to the nuclear or similar cases is not 
straightforward because i) the self-consistent mean field is absent, ii) the total energy is not conserved. 

A step forward in this direction has been performed by Bauer et al. in ref. [l5| ]. In this work the authors study 
the damping of collective motion in nuclei within the semiclassical Vlasov equation. Here selfconsistency is taken into 
account and the total energy is conserved. A multipole-multipole interaction of the Bohr-Mottelson type is adopted 
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for quadrupole and octupole deformation. In both cases the dynamical evolution shows a regular undamped collective 
motion which coexists with a weakly chaotic single-particle dynamics. 

In ref. [ fl6| we introduced a simple model which can shed light on the correlation between chaoticity and dissipation. 
We considered a gas of classical non-interacting particles moving in a two-dimensional billiard with nuclear dimensions. 
Particles collide with the oscillating walls and transfer energy to it, but the container can give back this energy heating 
the gas. We considered the gas + billiard as an Hamiltonian system, therefore the total energy is conserved with a 
good accuracy. Though no explicit dissipative term is considered the main effect of the coupling of the wall with the 
particles is a damping of the collective motion. The particle - wall coupling considered in our model is relevant for 
at least two other reasons : i) the coupling can enhance collectivity of the motion, since the particles are indirectly 
coupled among each others through the wall, ii) only with coupling included the motion of the particles can be driven 
to equilibrium at large timescales. In fact, it has been shown [ |l9[ ] that without the coupling the asymptotic particle 
velocity distribution is non-maxwellian. 

In ref. |l6[ ] only the monopole case was considered together with a small number of particles, namely 1 and 10. 
Here we extend the previous study by considering more particles (N=30) and L—2,3 multipolarities. For a fixed value 
of the wall deformation, the motion of the particles is regular in the case of the monopole L=0, partly chaotic for 
large deformation in the quadrupole case L=2, and essentially chaotic in the octupole case L=3. In this way we are 
exploring a set of physical conditions for which the motion of the fast degrees of freedom ranges from nearly ergodic 
to integrable (for a fixed value of the slow variable). The parameters of the model are chosen as typical at the nuclear 
scales. However, the considered set of dynamical systems should be generic enough to be representative of the physical 
problem under study, namely the damping process of a slow degrees of freedom in a bath of many fast degrees of 
freedom, far enough from adiabaticity and for different degrees of ergodicity. We therefore expect that the results we 
have obtained are generic enough to be qualitatively valid for a wide class of physical systems. 

The most important findings of our study are: a) chaos shows up in the single particle motion for any surface 
deformation for long time scales; b) the different geometry influences only the timescale for the onset of chaoticity 
which is faster for the higher multipolarities; c) the different timescales for the onset of chaoticity are in any case 
equal within a factor of two; d) the dissipation of the collective motion is sensitive only to this timescale and thus 
it depends only slightly on the geometry of the billiard for L=2 and L=3; e) though no explicit dissipative term is 
considered in our model, the damping of the collective motion is in practice irreversible due to the large Poincare' 
time. 

The present paper is organized as follows. We summarize the details of the model in section 1. The numerical 
results are illustrated in section 2. Conclusions are drawn in section 3. 



I. THE MODEL 

In ref. |ll| we considered a classical version of the vibrating potential model (see e.g. ref. [p0[ ). In this model several 
non-interacting classical particles move in a two-dimensional deep potential well and hit the oscillating surface. Using 
polar coordinates, the Hamiltonian depends on a set of {r^, 0j} variables, describing the motion of the particles, and 
the collective coordinate a. The Hamiltonian reads 

^22 2 

H(r iM = £(£ + + V(n, m) ) + fa + -MOV (1) 
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where {PniPOuPa} are the conjugate momenta of {ri,6i,a}. m = 938 MeV is the nucleon mass, and M = r/mNR^ 
is the Inglis mass, chosen proportional to the total number of particles N and to the square of the circular billiard 
radius R . The value of the factor rj is fixed in such a way to minimize the equilibrium fluctuations. In our case r] = 1 
for the monopole oscillation, whereas rj = 10 for quadrupole and octupole. Therefore in the L = case, collisions of 
particles against the walls are more inelastic, fi is the oscillation frequency of the collective variable a. The potential 
V(r,R(8)) is zero inside the billiard and a very steeply rising function on the surface, V(r,R(9)) = 



with V — 1500 MeV and a — 0.01 fm. Such a small value of the diffuseness was chosen in order to simulate 
closely a billiard system. Larger values of a should not affect qualitatively the results. The surface is described by 
R(0) = i? (l + a cos(LO)). Therefore this potential couples the collective variable motion to the particles' dynamics 
and prevents particles from escaping. The numerical simulation is based on the Hamilton's equations 

h = — , p ri = -^3 - Tp- (2) 

m mr^ or; 
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Please note that dropping the V term in the equation for p a , the wall motion becomes purely harmonic and decoupled 
by the particles' motion, as assumed in the model proposed by the authors of ref. |ll| , ^0[ . 

We solve the Hamilton's equations with an algorithm of fourth-order Runge-Kutta type with typical time step sizes 
of 1 fm/c. If not otherwise stated, the calculations were performed with the number of particles N = 30. The total 
energy was conserved with relative error AE/E < 10 -4 . For runs with the longest time duration or with the largest 
number of particles a fourth order symplectic integrator pl| ] was used, since it turns out to be more efficient with 
comparable degree of accuracy. 

As far as the initial conditions are concerned, we assign random positions to the particles inside the billiard 
and random initial momenta according to a two-dimensional Maxwell-Boltzmann distribution with a temperature 
T = 36 MeV. This value of temperature is chosen in order to mimic the Fermi motion of the particles. In this way 
the average velocity of the particles is close to the typical nuclear Fermi velocity. In a classical description the use of 
a Maxwellian distribution ensures that the initial conditions are not too far from equilibrium for the whole system. 

We consider the wall oscillations taking place not too far from adiabaticity. For this purpose we follow the definition 
of the adiabaticity parameter k given by the authors of ref. [fll|| , i.e. 

k = (5) 



being v the most probable particle speed, v — y/T/m, a a the initial amplitude of the oscillation and R the radius 
of the circular billiard. We choose R = 6 fm and O = 0.05 s _1 so that 

k = 1.53a G (6) 

With our choice of a (see below) , the adiabaticity parameter k might be not as small as assumed in the asymptotic 
theory of ref. The adopted value can correspond more to realistic situations. 

Since in the realistic cases the collective motion takes place around equilibrium, the initial wall coordinate has 
been chosen equal to ao = a + 6a, where a is the equilibrium value and 6a the deviation. The equilibrium value a 
corresponds of course to the thermodynamic limit, which is actually reached when considering an ensemble of copies of 
the system, all of them with an initial value a = a and differing from each other in the initial microscopic distribution 
of particles' positions and momenta. The equilibrium value depends on the considered multipolarity. In all cases we 
checked numerically that starting from the thermodynamical value a = a and a Maxwell-Boltzmann distribution, 
the collective variable a oscillates in time around a. Furthermore, we considered an ensemble of this type of runs, 
each one corresponding to different particles' initial conditions, consistent with a space uniform distribution and a 
Maxwell-Boltzmann velocity distribution. We then checked that the average over the runs of the value of collective 
variable a(t), at any given time t, was converging indeed to the constant equilibrium value a for a reasonable large 
number of runs. In this way one can also estimate the fluctuations around equilibrium the variable a undergoes in 
a typical run. Also the average amplitude of the fluctuations turns out to be consistent with the thermodynamical 
estimate for a harmonic oscillator in a thermal bath. More details on the procedure can be found in ref. (l6). After we 
checked that the numerical simulation produces good equilibrium properties, we perturbed the equilibrium collective 
coordinate a by an amount 6a = 0.3, and let both the billiard and the particles evolve in time. The chosen value is 
larger by about a factor of 3 than the average equilibrium fluctuations. Moreover at time t = we put p a = 0, the 
wall having only potential energy. 



II. NUMERICAL RESULTS 



A. Scatter plots 

One possible way in order to investigate the role played by the coupling and see whether it can induce chaotic 
dynamics, is drawing Poincare's surface of sections for the single particle coordinate. However this is impossible to 
perform in our case because of the large number of degrees of freedom. Then an alternative way to visualize a chaotic 
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behavior is to draw scatter plots, see Fig.l. There we display the final radial coordinate at a time t of one particle vs. 
the one at t = for the monopole, quadrupole and octupole deformation. The chosen times are t = 0.4r, 1.6r, 2.4r 
and 4r, r being the period of the oscillation, r = These plots are very useful and are commonly used in transient 



irregular situations like chaotic scattering |22f| and nuclear multifragmentation [|23| . The idea is that if the dynamics 
is regular, two initially close points in space stay close even at later times, but if the dynamics is chaotic the two 
points will soon separate due to the exponential divergence induced by chaos. In the first case this plot will show a 
regular curve, whereas in the other one a diffused pattern appears. We note that for all multipolarities the initially 
regular curves change into scattered dots, which clearly show that the coupling to the wall oscillation randomizes the 
single particle motion. This is at variance with what was discussed in ref. jllj], where chaos is supposed to appear 
only for multipolarities L > 2. In our model the coupling between wall and particles' motion produces a chaotic 
dynamics even for L < 2. In addition, the higher the multipolarity the earlier chaos starts because of the increased 
shape irregularity. 

Fig. i 
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FIG. 1. The final radial coordinate for one particle is drawn as a function of the initial one at different times t = 0.4t, 1.6t, 2.4t and 
At. Calculations are performed for multipolarities L = 0, 2, 3. 1000 initial conditions are considered. 
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Fig. 2 




FIG. 2. The time evolution of the fractal correlation dimension D2 is plotted for the multipolarities L = 0, 2, 3. The lines are to guide 
the eye. The dashed line is the fully random limit for scattered points on a plane. 

A more quantitative analysis can be performed because scatter plots of Fig.l have a typical fractal structure. As 
already done in rcf. [25(] , a fractal correlation dimension D2 can be calculated from the correlation integral C(r) 
One first counts how many points have a smaller distance than some given distance r. As r varies, so does the 
correlation integral C(r), defined as 

1 M 

where Q is the Heaviside step function and Zi a vector whose two components (xi,yi) are the points coordinates. M 
is the total number of points. The fractal correlation dimension D2 is then defined by 

InCtr) 

D 2 = lim — — — . (8) 

r— > In r 

Therefore by plotting the logarithms it is possible to extract D2 by fitting the linear slope for sufficiently small r. 
We considered as a good interval satisfying eq.(||) the one between r min = 1.83 • 1CT 2 and r rnax = 1. An ensemble of 
1000 points was considered. In Fig. 2 we display D2 vs. time for each multipolarity. At the very beginning D2 is equal 
to one, showing that the motion is regular. As time goes on, regularity is lost and the motion becomes chaotic until 
a complete randomness is reached, in which case D 2 = 2 as expected for a completely random distribution on the 
plane [ p3| . This result confirms the one published in ref. ]l6| ], where a different method of calculation was however 
employed. 



B. Lyapunov exponents 

A more common quantitative way to characterize the dynamics in a chaotic regime is by means of the largest 
Lyapunov exponent Ai, which gives the average rate of exponential divergence of two trajectories with nearly identical 
initial conditions. In general, if we denote with d(t) the distance between two phase space trajectories and d a = d(0) 
, Ai characterizes the average growth of the distance d(t) with time 

< d(t) > = d exp(Xit). (9) 
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Fig. 3 
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FIG. 3. For the cases L=0,2,3 we show the frequency distributions of the largest Lyapunov exponent Ai for 1 (thin histogram) and 30 
particles (thick histogram). An ensemble of about 40 events was considered in each case. 

We have calculated this largest Lyapunov exponent Ai by the method of Bcncttin et al. J2^] . This method consists 
in evolving two close trajectories originally separated in phase space by d 0l for a given time interval r, after which the 
magnitude of their separation d(r) is rescaled back to d a . The procedure is then repeated k times. It can be shown 
that Ai is given by the following limiting procedure : 

Ai = lim lim - — In— 7 — (10) 
fe^oo d o ^0 kr i d a 
i—l 

The result turns out to be essentially independent of d D . We considered the following metric in phase space 
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J2(s^ + s P l 3 + se] + spl) (11) 
3=1 



where the infinitesimal distances where normalized to the mean values and therefore are dimcnsionlcss quantities. 
In our case the length of the time interval r was 200 fm/c and d Q ~ 2.10~ 6 . Using the above method, a large 
number of trajectories (of the order of 40) was sampled for each multipolarity. Each trajectory was followed for a 
time t = 3500 fm/c, thus allowing for a safe determination of the largest Lyapunov exponent. 

In Fig. 3 we display for each multipolarity a distribution of Ai for a number of trajectories and for 1 and 30 particles 
(respectively thin and thick histogram). The most probable values do not depend strongly on the multipolarity, at 



variance with the results found in rcf. |11| and lie within a factor of 2 one from each other, the case L = being 
intermediate between L = 2 and L = 3. This is probably due to the different choice of the wall mass. In the case 
N=l, Ai is always smaller than the value obtained with N=30, thus showing that a) chaos starts earlier when a higher 
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number of particles is present and b) for long timescales the degree of chaoticity is more dependent on the number of 
particles than on the multipolarity. Therefore an Hamiltonian with coupling terms gives rise to a dynamics weakly 
dependent on the geometry of the billiard. It has to be noticed that the average number of particle-wall collisions per 
oscillation period is only about 2 for N=l and about 60 for N=30. 



C. Chaos and dissipation 



1. Analysis of the single event 



In this section we will mainly discuss the interplay between chaoticity of the single particle dynamics and dissipation 
of the collective motion. For this purpose, we follow at the same time the dynamics of the wall and the particles. 

Fig. 4 
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FIG. 4. On the left-hand side the time evolution of the collective variable is shown for the multipolarities L = 0, 2, 3 and on the 
right-hand side the corresponding excitation energy of the gas, normalized to the initial energy E . 

First we analyze the behaviour of one single event, keeping in mind that a correct statistical description can be 
performed only for an ensemble of events, as already pointed out in M]. On the left-hand side of Fig. 4 we plot 
the evolution of the square of the collective variable a 2 vs. time for one single event. Each panel corresponds to a 
fixed multipolarity. We note that the amplitude of the collective motion shows an irregular oscillation, at variance 
with the results found in rcf. |l5|. We note also that a slight damping can be observed, in contrast with the results 
published in ref. |Iq| , where calculations with a smaller number of particles were performed. Please note that a 2 
keeps on oscillating around the equilibrium value a which for the monopole case depends on the gas temperature and 
on the wall frequency |{iq| . For the other multipolarities a oscillates around zero and a 2 decreases more rapidly for 
increasing L, indicating a faster dissipation for a larger irregularity of the billiard. 
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Fig. 5 
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FIG. 5. The square of the ratio a/a is plotted vs. the number of oscillations for different multipolarities and number of particles. 



On the right-hand side of Fig. 4 we plot the time evolution of the excitation energy of the gas, defined as the relative 
variation of the total energy of the gas E with respect to its initial value Eq, AE/Eq. In all three cases the gas is 
heated up, but with a different time dependence according to the multipolarity. Except for small irregular fluctuations, 
an increasing trend shows up for the quadrupole and octupole modes. An oscillating pattern, slightly irregular, is 
visible for the monopole case. Moreover, in the L = oscillation the energy gained by the gas is lower than in the 
L = 2, 3 cases. Therefore some dissipation is present for all multipolarities and is larger for increasing L. Of course 
these are only general features of the events for different multipolarities, the details being different from one event to 
the other. 

In order to better understand the mechanism of dissipation, we display in Fig. 5 the square of the ratio a/a , being 
a D = a(t = 0). The plot is as function of time, for different multipolarities L and particle number N. In our model 
the total energy is conserved so the damping of the wall motion corresponds to the heating of the gas. This process is 
in principle reversible, but the high number of degrees of freedom involved makes in practice this process irreversible. 
One should wait a very long time, i.e. the Poincare' time, to see the system again close to the initial conditions. We 
note that the dissipation is stronger when both L and N increase. In particular, the increased number of particles 
makes the available phase space bigger and the Poincare' time longer. The latter can therefore be a limitation, in the 
sense that a short Poincare' time would weaken dissipation. 



2. Analysis of an ensemble of events 



In order to have a global picture of the macroscopic system, many events are needed and average values of the 
different obscrvables should be considered. In Fig. 6 we display the time evolution of a 2 averaged over an ensemble 
of 1000 different events, each obtained by assigning random initial conditions to the particles both in coordinate and 
momentum space but consistent with a Maxwellian. We note that the motion of the collective variable < a 2 > is 
completely damped out for times which depend on the multipolarity and develops around equilibrium, as we checked 
explicitly. The reader should keep in mind that for L = collisions of the particles with the wall are more inelastic (see 
paragraph 2) due to the lighter mass of the wall, therefore the monopole oscillation damps out in a time comparable 
with the one of the L = 2 mode. The same mass was used in the calculations published in ref. If we put 77 = 10 
even for L — 0, the damping time is longer. 
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Fig. 6 




FIG. 6. The square of the collective variable a, averaged over an ensemble of 1000 events, is displayed vs. the number of oscillations 
for the three multipolarities considered. 
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FIG. 7. On the left-hand side the wall energy, averaged over an ensemble of 1000 events and normalized by its initial value, is reported 
vs. the number of oscillations for different multipolarities. On the right-hand side the observable Ao 2 is drawn (see text for further 
details). 

The damping of the oscillations, which is apparent in Fig. 6 for all the three multipolarities, can have two distinct 
origins. On one hand, each individual event can indeed display a damping of the oscillation amplitude. In this case 
the effect of the average is only of smoothing the fluctuations. The observed damping time is then simply the average 
damping time of a generic set of events. On the other hand, each event has a different time evolution, and after a 



9 



period of time the phase of the oscillation can be quite different from one event to another. If the phase becomes 
essentially random as the time proceeds, when taking the average of a over different events strong cancellations can 
occur and a damping of the oscillations can appear, even in the case where no average energy damping is present. 
It is worth noticing that this damping mechanism disappears if the ratio between the single particle mass and the 
wall mass is vanishing small, and therefore it is ineffective at macroscopic level, where event-to-event fluctuations can 
be neglected. In order to estimate the contribution of each one of the two mechanisms to the damping displayed in 
Fig. 6, we plot on the left-hand side of Fig. 7 the energy of the wall divided by its initial value, averaged over the 
same sets of events. As a function of time, this quantity should be insensitive to phase randomization because it is 
proportional to < a 2 > and corresponds to the average energy damping of the oscillations. 

For comparison, on the right-hand side we plot also Aa 2 = (< a > -doo) 2 , a quantity which, on the contrary, 
should be sensitive to phase cancellations. For the monopole case, the decay time of the wall energy to its asymptotic 
value is substantially longer than the one for Aa 2 , which indicates that for the monopole case the dephasing 
mechanism is quite important. The asymptotic value of < a > is consistent with the equilibrium value a for a gas 
at the final temperature, which, from energy conservation, turns out to be T = 43 MeV. We checked that indeed 
the velocity distribution, averaged over the events, is very close to a Maxwellian with that temperature. Furthermore 
the average amplitude of the fluctuations around the equilibrium value was checked to be consistent with a harmonic 
oscillator in equilibrium with a thermal bath at T = 43 MeV. 

On the contrary, for the higher multipolarities the decay time of the wall energy and Aa 2 appear to be very close, 
and therefore the dephasing mechanism seems to be in these cases ineffective. The damping is mainly an energy 
damping and is characteristic of a generic event. 

Fig. 8 
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FIG. 8. The excitation energy of the gas, averaged over an ensemble of 1000 events, is reported vs. the number of oscillations for 
different multipolarities. 

These considerations are confirmed if one calculates the excitation energy of the gas, displayed in Fig. 8. We note 
that the characteristic times for the energy growing follow closely in all cases the decay times of the wall energy. 
Moreover for L = 2,3 we note the presence of two different regimes: a first one lasting for about one two oscillations 
characterized by a sharp rising in the excitation energy, and a second one at successive times where a slower approach 
to equilibrium is apparent. As it can be deduced from Fig.l, the first stage may have some relation with the onset of 
chaos in the single-particle motion. After chaos has fully developed, and an appreciable part of the total energy has 
been pumped into the gas, a certain degree of randomization is reached and a slower dissipation rate shows up. A 
similar behavior in the relaxation was observed for the HMF model [ p~7j and the hard sphere gas and the Lorentz gas 
studied in ref. p8[ . In the latter cases the time evolution of the Boltzmann entropy was studied. 

This behavior is completely absent in the monopole case. Therefore it seems that different kinds of dissipation can 
originate from the same underlying chaotic single-particle motion. It should be stressed that, within the time of chaos 
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development, the wall has dissipated only a fraction of its energy, as it is clearly seen in Fig. 7. In the monopole case 
it is apparent that the flow of energy from the collective variable a to the particle gas is not completely "irreversible" , 
and the wall receives back part of its energy from the gas for a number of oscillations. 

3. Comparison with the wall formula 

Let us now try to make some connection to one body dissipation and the wall formula introduced and studied 
extensively in the past by Swiatecki and coworkers jn]|nj. We limit ourselves to the analysis of the octupole mode, 
because for that (and for higher multipolarities) the Swiatecki model predicts dissipation. Of course a comparison 
between the situation studied here and in JTl| may be biased because in the model studied in there is no energy 
conservation. Therefore phase space of the gas is conserved in the latter model but not in the case considered here. 
In spite of that our model should be able to recover the Swiatecki limit. And actually it does. This happens when in 
our model the mass of the wall becomes very large as for example in Fig. 9 where we choose r\ = 100 and a = 0.09 
with N = 30 and £1 = 0.05 as before. In this case the particle-wall collisions are more elastic and the time evolution 
is slower for all observables. In particular in panels (c,d) of Fig. 9 the dots are our numerical simulation and the solid 
and dashed lines are two versions of the wall formula to be explained in the appendix. In particular, the dashed line 
represents the results of the wall formula which contains only terms linear in time. The inclusion of terms quadratic 
in time, which take into account for the increased average speed of the gas for large times, produces a new dissipation 
formula represented in the figures below by a solid line. Now, if every thing is consistent we should get the same 
result as the one shown in Fig. 9 in switching off the coupling. That this is actually the case is shown in Fig. 10 where 
the time evolution of the different quantities is practically the same as in Fig. 9. 



Fig. 9 



1.5 I ■ ■ ■ ■ i ■ ■ ■ ■ I 1.5 




T 



0.4 - (c) (d) 




: , , , , i , , , , i t , , , , i , , , , : 

2 4 1 2 

t/r t/r 



FIG. 9. Coupled octupole motion of the wall and gas particles with a very large inertia of the wall (rj = 100). The adiabaticity 
parameter is k = 0.14. Shown are the elongation (panel (a)), the wall energy (panel (b)), the relative increase of the gas energy (panel(c)), 
and a blow up of panel (c) (pancl(d)). 
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Fig. 10 




FIG. 10. Same as Fig. 9 but without coupling. The frequency of the wall is kept fixed to the initial one, i.e. f! = 0.05s 1 . 

A particular feature of the case shown in Figs. 9-10 is the fact that the wall formula results follow the exact 
evolution only during the first period. We should emphasize that we applied here the wall formula adapted to the 
two dimensional case, as derived in the appendix. In Fig. 11 we show a repetition of the calculation in ref. [p"i"| where 
the wall formula follows the exact result for the heating of the particles gas over several periods. The only difference 
with the calculation in [ [llfl is that we here use again for the gas particles our Maxwell-Boltzmann distribution with 
T = 36 M eV, whereas in ]ll| a sharp Fermi-Dirac distribution was employed. The fact that the wall formula for this 
particular case of the parameters is well reproduced indicates that this feature is independent of Maxwell-Boltzmann 
or Fermi-Dirac statistics. This must be expected, since in the wall formula the dissipation rate depends only on the 
average velocity of the particles. Therefore, the dissipation must be the same as long as the average velocity of the 
particles is the same. The only variance with the case shown in Figs. 9-10 is that the parameters are slightly different, 
or that for instance the adiabaticity parameter k defined in eq.(6) is by almost a factor of 4 smaller in Fig. 11 than 
in Figs. 9-10, that is k — 0.04 instead of re = 0.14. However, both values of k are sufficiently small in order to verify 
the adiabaticity criterion (for the validity of the wall formula) established by Swiatecki et al., namely that k << 1. 
According to our findings here, re must really be extremely small so that the wall formula holds during several periods, 
otherwise it may be valid only during the first instances of the dynamics. 
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Fig. 11 




FIG. 11. Repetition of Swiatecki et al. calculation but with a Maxwellian distribution for the gas particles. The observables displayed 
are the same as in Figs. 9-10. The adiabaticity parameter is k = 0.04. See the text for more details. 

Let us now come back to a more detailed analysis of our present results obtained including the coupling terms in 
the Hamiltonian. In Fig. 12 we show as a function of time the evolution of the elongation (panel (a)), the wall energy 
(panel (b)), the gain in energy of the gas (panel (c)) together with a blowup of (c) (panel (d)). In panels (c) and (d) 
are also shown the results of the standard wall formula (broken line) , which is linear in time |1 1| , and the new version 
including a quadratic term (full line) p2| . We should stress that we applied here generalized wall formulas taking 
account of the fact that the wall motion is actually damped. This implies that we have to use a time dependent 
version of the wall formula (TDWF) as derived in the appendix. 



Fig. 12 
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FIG. 12. Same as Fig. 9 with coupling. In panels (c) and (d) time dependent wall formulas have been used, taking into account 
respectively terms linear in time (broken line) with quadratic ones (full line). The dots are our numerical results. 

Again we see that the wall formula only agrees with the exact solution during the hrst period. For longer times the 
wall formulas give rise to much more energy dissipation than in the exact evolution, the new formula more than the 
standard one. 

From the above investigations it seems therefore that the wall formula can only simulate the dissipation of energy 
in the first instances of the dynamics when the adiabaticity parameter k (eq.(6)) is extremely small (of the order of 
1/100). 

III. CONCLUSIONS 

In conclusions, we have presented a dynamical approach based on the solution of the Hamilton's equations for 
several particles moving in a classical billiard having nuclear-like dimensions, in order to explain dissipation of the 
collective motion. We found that the presence of a coupling term in the single particle Hamiltonian induces chaotic 
motion at the microscopic level. 

As far as the monopole mode is concerned, we found irregular behavior together with damping in single events. 
The damping observed when an average is taken over a set of events is a consequence also of the irregular time 
dependence of the oscillations. This incoherence is produced by the chaotic single particle dynamics, which makes all 
events belonging to the same ensemble different one from each other, and therefore substantial cancellations occur 
once the average is taken. For the higher multipolarities this mechanism seems to be ineffective, and the damping 
coincides with a real energy damping of the oscillations due to the large Poincare' time. The dissipative process for 
the quadrupole and octupole modes looks different also in another respect. In fact, while the single event properties 
are similar to the monopole case, an ensemble of events shows that two different regimes appear : a) an initial fast 
dissipative evolution corresponding to the onset of chaos in the single-particle motion and b) a slower dissipative trend 
towards equilibrium. 

All these results should be qualitatively independent on the particular nuclear dimensions we have used in the 
model. 

In order to be closer to equilibrium we have used for the initialisation of the gas particles a Maxwell-Boltzmann 
distribution with a temperature of T = 36 MeV to have a mean kinetic energy characteristic of nuclear systems. 
We also checked that the use of a Maxwell-Boltzmann distribution versus a Fermi-Dirac step is apparently of no 
consequence on the damping rate. Indeed in switching off the particle wall coupling and periodically octupole- 
deforming the wall under the same condition as was done some time ago by Swiatecki and collaborators [jlTj| using 
a Fermi-Dirac step, we find with a Maxwell-Boltzmann distribution an identical rate of feeding of energy into the 
particle gas. Moreover we also checked that the wall formula prediction for energy dissipation agrees with the numerical 
simulation over several periods in time. This shows that our model is able to mock up real nuclear dynamics together 
with its damping mechanisms. 

When we increase the adiabaticity parameter k by almost a factor of four attaining k — 0.14, what should still 
be considerable as small, the agreement of the numerical result with the wall formula was reduced to only the first 
period of the oscillation. In view of what we said above about Maxwell-Boltzmann versus Fermi-Dirac distributions, 
we assume that this is a generic result. Also when reducing the inertia of the wall by a factor of ten and thus 
increasing the particle wall coupling when reestablished the heating rate of the particle gas agrees with the wall 
formula only during the first period. This happens in spite of the fact that we used a time dependent version of the 
wall formula where the continuous slowing down of the wall motion has been taken into account. Indeed the wall 
formula considerably overestimates the damping rate at longer times. 

One should realise that our model implies real particle- wall collisions with global energy conservation. Such collisions 
are absent in pure Hartree-Fock or Vlasov calculations. In fact our model in what concerns the particle wall collisions 
should come quite close to the situation considered e.g. in [3j where the collision term is based on the particle vibration 
coupling model. Also in that model the damping rate strongly differs between low and high multipolarities of the 
collective motion. However in ref. ||, contrary to the spirit of the present work, fast (diabatic) motion (the giant 
resonances) was considered, what makes a detailed comparison inadequate. In the future, in order to simulate a pure 
mean field dynamics based on our model, we plan to calculate the evolution with the so called parallel ensemble 
technique and ensemble average at each instant of time the motion of the wall. We anticipate that this very much will 
suppress the particle wall damping mechanism so that the collective oscillation will be mostly undamped, a feature 
which should bring the present study in closer contact with the one performed in ref. | fl5|| . 

Our model certainly has implications beyond nuclear physics, it should be generic for all situations where a heavy 
particle (here the wall) is moving in a Knudsen gas of light particles. For example the motion of a pendulum of mass 
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M in a very rarefied gas of particles with mass m << M should show similar features as the ones studied within the 
present model. If the big particle can be approximated by a sphere suspended on a spring and the gas is enclosed in 
a box, the whole system can be considered as some sort of dynamic generalisation of a Sinai billiard. The study of 
the damping of the heavy particle in such a situation would be particularly interesting. 
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APPENDIX A: 



Here we give some details about the two-dimensional wall formula. Starting from the original formulation given by 
Swiatecky in ref. [flo|| , after rescaling in two dimensions, we get at the following result 

dE 4 / 

— — = — pv i> q 2 dl + ....{terms of higher order in q 3 ) (Al) 
dt n J 

p and v are respectively the density and the average speed of the gas, q is the speed of the wall and dl is the line 
element. By integration of eq.(Al), and neglecting corrections of order a 2 or higher, we easily get 

^=4pvR 3 (a(t)) 2 (A2) 

where a(t) can be calculated either without or with coupling. In the former case a(t) = a cos(flt) and Eq.(A2) will 
read 

1 dE = 2Ro n 2 a 2 sln 2 m (A3) 

E dt v 

E and a are the initial energy of the gas and the initial elongation. After more algebra we get the excitation energy 
of the gas AE/E Q 

AE R n 2 r ^ sm(2M. ,. . . lS 

= 11 a a [lit no coupling (A4) 

E Q v 2 

In the case with coupling a(t) comes directly from our numerical simulation. In fact, since the wall energy is 

E w = T^Mct 2 + ^MVL 2 a 2 (A5) 
we easily get a and substitute in Eq.(A2) obtaining 

= — E w - SI a with coupling (A6) 

E dt v M 

By integration on time we get again the excitation energy of the gas AE/E a . The last equation is the time 
dependent wall formula (TDWF). 

In ref. ]l2[| , a generalized wall formula has been studied, which takes into account for the increased average speed of 
the gas for large times. The new dissipation formula contains an additional term quadratic in time, besides the linear 
term (eq.A4). We have performed similar calculations in two dimensions, both without and with coupling. The new 
formula for the excitation energy of the gas reads 

^=21{t) + ~m (A7) 

E 7T 



I(t) being respectively 



„ a 2 flR .^ sin(2Qt), , 

lit) = — — -Kit K - -} no coupling A8 

2v 2 



15 



R f 2 

I(t) = -1 \ dt[—E w - n 2 a 2 } with coupling (A9) 
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